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GLOBAL WELL-POSEDNESS OF A NON-LOCAL BURGERS 
EQUATION: THE PERIODIC CASE 

CYRIL IMBERT, ROMAN SHVYDKOY, AND FRANCOIS VIGNERON 


Abstract. This paper is concerned with the study of a non-local Burgers equation for 
positive bounded periodic initial data. The equation reads 

Wt-u|V|M+|V|(w^) =0. 

We construct global classical solutions starting from smooth positive data, and global weak 
solutions starting from data in L°°. We show that any weak solution is instantaneously 
regularized into . We also describe the long-time behavior of all solutions. Our meth¬ 
ods follow several recent advances in the regularity theory of parabolic integro-differential 
equations. 


1. Introduction 

Dynamics of fluid motion provide a rich source of evolution laws that defy a complete 
well-posedness theory. Apart from the classical three dimensional Euler and Navier-Stokes 
equations, scalar models such as the super-critical surface quasi-geostrophic (SQG) equation 
or Darcy’s law of porous media pose core difficulties, both to the traditional and newly 
developed approaches. In recent years, several classes of models have appeared that either 
mimic the basic structure of the aforementioned equations or serve as viable models on their 
own. For example, in relation to the SQG model 

et + u-ve = o 

where the div-free velocity u = R^9 is the perpendicular Riesz transform of 6*, A. Gordoba, D. 
Gordoba, and M. Fontelos [8] have studied an analogous ID model given by 6t+ {H6)6x = 0, 
where H denotes the Hilbert transform. SQG written in divergence form 

Ot + div(M0) = 0 

received its counterpart in the form of 9^ + (9H9)x = 0 which was studied in [5], as well as 
convex combinations of the two ID models above. The latter combination Erst appeared in 
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physics literature [1] as a model for a conservation law such as the classical Burgers equation 

lit ~ 

but with a non-local flux F = OHO. In these models the relation between the drift and the 
driven scalar is ensured by a Fourier multiplier with an odd symbol. This sets them apart 
from other laws such as porous media, Burgers or Euler equation or the recent Moffat’s 
scalar model of magnetostrophic turbulence of the Earth’s fluid core [22], In all the later 
cases, such a relation is furnished via an even symbol which generally entails a more singular 
behavior. 

In this paper we study the following model: 

(1) Mt - n|v|M + |V|(m^) = 0, a; € or T'’*, 

where |v| = denotes the square root of the Laplacian and has the symbol |.^|. With the 

addition of viscosity, the model ut + u\\j\u-\v\{v?) = i^Au was proposed by P.G. Lemarie as a 
scalar case study of the 3D Navier-Stokes (note the opposite signs). The works of F. Lelievre 
[18], [19], [17] presented the construction of global Kato-type mild solutions for initial data 
in and of global weak Leray-Hopf-type solutions for initial data in 

and A local energy inequality obtained for this model was suggestive of possible 

uniqueness for small initial data in critical spaces, in a similar fashion to 3D Navier-Stokes. 
The focus of the present paper will be on the inviscid case z/ = 0. 

Without viscosity the model bears resemblance to some of the “even” inviscid cases in the 
sense explained above. For example, dropping the 1/2 factor, the Burgers equation can be 
written in the form of a commutator: 


dtu = [u,dx]u. 

Our model replaces with the non-local operator |V| of the same order (hence our choice 
of name for (1)). The classical incompressible Euler equation is given by 

Ut + u- Vu + Vp = 0, 

where p is the associated pressure given by p = T{u0u) + local, where T is a singular integral 
operator with an even symbol. We thus can draw an analogy between terms: u-Vu —ulvlu 
and Vp ~ |v|(m^). Actually, analogies with Euler or Burgers extend beyond the formal 
range. Let us discuss the basic structure of (1) in greater detail. In what follows we give our 
model (1) the name of the Non-local Burgers equation. 

At a formal level, (1) shares a more intimate connection with other equations of fluid 
mechanics. For example, if one applies formally the commutator theorem in the scalar case, 
even though u might not be smooth and |.^| surely isn’t, one gets: 

[u, |V|] = Op u'(x) + 5^=0 u"(x) + ... j . 

Thus in dimension 1, the equation becomes formally dtU = u'Ru + {J^u)u" + ... where R 
denotes the Riesz transform. The term u'Ru is a scalar flavor of the SQG nonlinearity 
written above. 
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Let US recall that, in the operator |v| enjoys an integral representation in terms of 
convolution with the kernel K{z) = for some normalizing constant Q > 0 depending 
only on the dimension. The model (1) can thus be rewritten in an integral form: 

(2) dtu = p.v. / K{y-x)(u{y)-u{x))u{y)dy. 

If u is periodic with period 27t in all coordinates, the representation above can alternatively 
be written as 

(3) dtu = p.v. / Kper{y-x){u{y)-u{x))u(y)dy, 

where is the torus and iLper(^) = More explicitly in ID, iLper(^) = 

For periodic solutions, both representations are valid due to a sufficient decay of K at inhnity. 
The former is more amenable to an analytical study due to the explicit nature of the kernel 
and applicability of known results, while the latter will be more useful in our ID numerical 
simulations presented at the end. 


The following basic structure properties of the model can be readily obtained from either 
representation, at a formal level. Let u be a solution to (1). 

(TI) Translation invariance: if Xq e to > 0 then u{x + Xo,t + to) is another solution. In 
particular, the periodicity of the initial condition is preserved. 

(TR) Time reversibility: if to > 0, then -u{x,to - t) is a solution too. 

(SI) Scaling invariance: for any A > 0 and a,/3 € M, X°‘u{X^x, X°‘'*'^t) is another solution. 
(MP) Max/Min principle: if u > 0, then its maximum is decreasing and its minimum is 
increasing. 

(AMP) Anti-Max/Min principle: if n < 0, then its maximum is increasing and its minimum 
is decreasing. 

(E) Energy conservation: ||M(t)||j ;,2 = ||mo||l 2 is obtained by testing (1) with u. 

(HP) Higher power law: for any p e (2, oo), the following quantity is conserved: 

\\u(t)\\l^ ^ 2 X Jln^ Hx)u(y){\u{y)\P-^ - \u{x)\P-^){u(y) - u{x))K{x,y)dxdy. 

This equation, obtained by testing (1) with \u\p~‘^u, implies the decay of those Lp 
norms when u > 0. 

(ML) First momentum law: for either D = or = T'^, integrating (2) and using the 
Gagliardo-Sobolevskii representation of the iL2(y2)-norm (see [9]): 

(4) / u(x,t')dx = / u(x,t)dx+ / ||m(s)P. 1 ds. 

Jn Jn Jt H^(n) 

On T'^, one has c and this property combines nicely with (E) and ensures that 


regardless of the sign of uq- 

For positive solutions u > 0 the right-hand side of (2) gains a non-local elliptic structure of 
order 1, while for u < 0 the equation behaves like a backward heat equation. However, with 
the energy conservation (E), the model shares common features with conservative systems 
such as Euler, giving it a second nature. As a result, for u > 0, we see an accumulation of 
energy in the large scales and a depletion of energy on small scales. More specihcally, as 
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seen from from the momentum law (ML), the dissipated energy from high frequencies gets 
sheltered in the first Fourier mode. At least on the phenomenological level, this property 
parallels what is known as the backward energy cascade in the Kraichnan theory of 2D 
turbulence (see [16] and references therein). This makes our model potentially viable in 
studying scaling laws for the energy spectrum, structure functions, etc. 


The aim of this paper is to develop a well-posedness theory for the model and study its 
long-time behavior. We make use of both sides of its dual nature by way of blending classical 
techniques relevant to the Euler equation, such as energy estimates, a Beale-Kato-Mayda 
criterion, etc [20], with recently developed tools of regularity theory for parabolic integro- 
differential equations [4, 7, 10, 12, 13, 21]. Let us give a brief summary of our results with 
short references to the methods used. Complete statements are given in Theorems 2.1, 2.2, 
2.5, 2.6, 3.2. First, we declare that all the results are proved in the periodic setting, except 
local existence which holds in both the periodic and the open case. Periodicity provides extra 
compactness of the underlying domain which for positive data, due to the minimum principle 
(MP), warrants uniform support from below in space and time, which further entails uniform 
ellipticity of the right-hand side of (2). 


Local existence. For initial data Uq e on or T'^, with mq > 0 pointwise and 

m > d/2 + 1, there exists a local solution in C'([0, T); iL'"(D'^)) n C^([0, T); i7"*“^(D'^)). Even 
for this local existence result, the positivity of the initial data is essential. We also have a 
BKM regularity criterion: if || Vn(t) ||L<>=dt < 0, the solution extends smoothly beyond T. 
The proof goes via a smoothing scheme based on a desingularization of the kernel. 

Instant regularization. Any positive classical solution to (1) on a time interval [0,T) 
satisfies the following bounds: for any k e N, there exists ak e (0,1) such that for any 
0 < to < 

||M||(^fe+-fc.-fe(Trdx(to,T)) - C{d,k,to,T,mmuo,m&xuo) 

(5) 

l|V^;«||Lip(Tdx(to,T)) ^ C{d,k,to,T,mmuo,m&xuo). 


To achieve this we symmetrize the right-hand side of (2) by multiplying it by u and then 
writing the evolution equation for w = u^\ 

(6) dtw = p.v. J {w{y) - w{x))k{t, x, y)dy 

The active kernel k is symmetric and satisfies uniform ellipticity bound < k{z) < 

This puts the model within the range of recent results of Kassmann et al. [2, 14] and of 
Caffarelli-Chan-Vasseur [4] where Moser / De Giorgi techniques were adopted to yield initial 
Holder regularity for w and hence for u by positivity, i.e. estimate (5) with the index /c = 0. 
See also [15, 6]. We then apply our new Schauder estimates for parabolic integro-differential 
equations with a general kernel [12], see also [13, 21], to obtain the full range of bounds (5). 
Global existence. It readily follows from the BKM criterion and the instant regularization 
property. 

Global existence of weak solutions. Since the bounds (5) depend essentially only on the 
L°°-norm of the initial condition, we can construct global weak solutions starting from any 
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Mo e Uq > 0. The solution belongs to the natural class T);nL°°([0,T) x 

T'^) for all T > 0. The initial data Mq is realized both in the sense of the L°° weak*-limit and 
in the strong topology of Such a solution satisfies (5) instantly. A surprising difficulty 
emerged here in recovering the initial data since the sequence of approximating solutions 
from mollified data may not be weakly equicontinuous near time t = 0. A weak formulation 
of (1) does not allow us to move the full derivative onto the test function. See Section 2.4 
for a complete discussion. 

Finite time blowup for uq < 0. As follows from the previous discussion and (TR), a 
negative solution may start from C°° data and develop into L°° in hnite time. We don’t know 
whether a more severe instantaneous blowup occurs for negative data. Our requirement for 
local existence certainly supports this. However, the numerics presented in Section 5 suggests 
that for some mostly positive data with a small subzero drop, the positive bulk of the solution 
may persevere. The solution gets dragged into positive territory and exists globally. 

Long-time asymptotics. As t ^ +oo, any weak solution to (1) converges to a constant, 
namely consistent with (E), in the following strong sense: the amplitude of u{t) tends 

to 0 exponentially fast with a rate proportional to minwo- The semi-norm ||VM(t)||Loo does 
the same. Thus, there are no small-scale structures left in the limit. The latter statement 
requires more technical proof which relies on an adaptation of the recent Constantin-Vicol 
proof of regularity for the critical SQG equation, [7]. 

The organization of the paper follows the order of the results listed above. To shorten the 
notations, we frequently use || • ||m to denote the Sobolev norm of and | • |p to denote the 
LP-noim. 


2. Global well-posedness with positive initial data 

2.1. Local well-posedness. We start our discussion with local well-posedness in regular 
classes. Let denote or T'^. 

Theorem 2.1 (Local well-posedness). Given a pointwise positive initial data uq € 
where m>(i/2 + l is an integer, there exists a time T > 0 and a unique solution to (1) with 
initial condition uq, which belongs to the class G([0,T); 77^(12'^)) n C'^([0,T);77™'“^(r2'^)). 
Moreover, u(x,t) > 0 for all {x,t) € x [0,T), and the maximum maxQdM(t) is strictly 
decreasing in time. 

The proof in the case of M'’* requires slightly more technical care in the maximum principle 
part, while being similar in the rest of the argument. We therefore present it in only. 
In the case of the torus however, we will also obtain a complementary statement for 
the minimum: minTrdM(t) is a strictly increasing function of time, thus the amplitude is 
shrinking. In Section 3 we will elaborate much more on the asymptotic behavior of the 
amplitude. 

Proof. The proof will be split into several steps. 

Step 1: Regularization. Let us consider the following regularization of the kernel 

Cd 
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and the corresponding operator 

|V| 5 M= / Ks{x-y){u{x)-u(y))dy = u Ks{z)dz-Tsu 
where Tg is the convolntion with Kg. The regnlarized eqnation takes the form 

(8) dtu = [u,\\/\g]u = / Kg(x-y){u{y)-u(x))u(y)dy = -[u,Tg]u. 

JR'* 

Note that Tg is inhnitely smoothing, i.e. < cys,s'||M||s', for any 0 < s' < s. So, by 

the standard qnadratic estimates, the right-hand side of (8) is qnadratically bonnded and 
locally Lipschitz in Thns, by the Fixed Point Theorem, there is a local solntion u e 
C^([0,T);iJ™') with the same initial condition uq. Here T depends on ||Mo||m and 6. For 
later use, note that \u(t )\2 = |mo |2 is conserved. 

Step 2: Maximum principle. Suppose Uq e iJ*", uq > 0, and u e C^([0, T); JJ™) is a local 
solution to (8) with the initial condition uq. As ^ for m > (i/2+1 (even better 

is true), and u(x,t) 0 as x oo, then u{t) has and attains its maximum M{t) = maxjjd u{t). 
We claim that u{x,t) > 0, for all {x,t) € x [0,T), and the maximum function M(t) is 
strictly decreasing on [0,T). Let us prove the hrst claim hrst. 

Let us £x i? > 0 and show that u never vanishes on (0,T) x Suppose it does. Let 

us consider 

to = inf{t 6 (0,T) : 3|x| < R,u{x,t) = 0}. 

By the boundedness of (0,T) x 5/^(0) and the continuity of u, R is attained. Thus, since 
Uq > 0, then to > 0. Let xq e 5^(0) be such that u{xo,to) = 0. Evaluating (8) at (xo,to) we 
obtain 

Ut(xo,to)= / ^s(xo-y)u^(y)dy>0. 

jRd 

Observe that the right-hand side is strictly positive since the energy of solutions to (8) is 
conserved. This shows that for some earlier time t < to there exists x e 5^(0) where u 
vanishes, which is a contradiction. Since the argument holds for all i? > 0, the claim follows. 

Let us prove the second claim now. Suppose that M{t) is not strictly decreasing on 
[0,T). This implies that there exists a pair of times 0 <t' <t" <T such that M{t') < M{t''). 
Let us show that there exists a to > t', such that M(to) > M(t) for all t e [t',to]- If 
M(t') < M{t”), then by the continuity of M{t), M attains its maximum on the interval 
[t',t"]. Let to e be the left outmost point where the maximum of M is attained. 

Then to > t', and M(to) > M(t) for all t' < t < to. If, on the contrary, M{t') = M{t") then 
either one can shrink the interval to fullhll the previous assumption or M{t) is constant 
throughout In either case, there exists, as claimed, ato> t', such that M(to) > M(t) 

for all t' <t < to. Let us consider a point xq € such that u(to,xo) = M(to). Then 

(9) dtu{to,xo)= Kg{x-y){u(y)-u{xo))u{y)dy <0. 

jRd- 

So, at an earlier time t < to, one must have u{t,xo) > M(to) in contradiction with the initial 
assumption. 

Step 3: 6-independent bounds. Let us observe the following representation formula that 
follows easily from u{x) - u{y) = Vm(( 1 - A)x + Xy) ■ {x - y)d\ : 

|V|5M = f Rx5{Vu)dX, 

Jo 
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where is the smoothed (vector) Riesz transform given by a convolution with the kernel 
Let us notice that $ 7 (^) = (-), while is the z-multiple of the 

(72+|2p)^2^ 7 V '/ 

Poisson kernel. Therefore, on the Fourier side, 


$7(0 


lel • 


Since these symbols are uniformly bounded, the family of operators is uniformly 

bounded in and any Sobolev space H^. Thus, we have uniform estimates 


(10) \\Msu\\k<C{d,k)\\Vu\\k, 

for all k >0. Finally, note that, in any space dimension, the full symbol of |v |5 is given by 

(11) Sym|v|5 = |e| re-^^l«ldA = ^(l-e-^l«l). 

Jo 0 


Let s be a multi-index of order |s| = m. Differentiating ( 8 ), we obtain 

dtd''u= Y, d"~°‘u\v\5d°‘u + u\v\sd'^u-2\v\5{ud"u)- Y |V|5(9“Ma^"“M). 


0<Q:<S 

Let us test with d'^u. We have: 
1 
2 


0<a<s 


( 12 ) 


dt\d^u\l= f d^u ^ °‘u\v\5d°'u+ f ud^u\v\5d^u-2 f d^u\v\diud^u) 

0<a<s d J 


- f d^u ^ | v | 5 ( 9 “ m 9 '^ “ m ). 

0<Q:<S 


The middle two terms of (12) contain derivatives of order m + 1. By symmetry however, 
they add up to -/ ud^u\\/\sd^u. By the positivity of m, a bound for this term follows from 
the elementary identity -a(a -b) < - 5 ( 0 ^ - b'^): 

- J ud^u\v\sd^u = - JJ u{x)d^u{x){d^u{x) - d^u{y))Ks{x - y)dxdy 

< -^ JJ u{x){d^u{xy - d‘^u(yy)Ks{x - y)dxdy 

= / uM,(a-uy = -1 J(M,u)(d‘uy. 

One thus gets, using (10) for the last step: 

- J ud'^ulvlsd'^u - ^ J (|V|5m)(9^m)^ < ||V|5M|oo||M||m ^ ^ ||M||m 

where all the constants in the inequalities are independent of 6. 

The rest of the expression (12) is simpler to deal with as it does not contain any other 
derivatives of order m +1. To estimate the hrst sum we use the Gagliardo-Nirenberg inequal¬ 
ities: 
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So, for any 0 < a < s, we have 


/ 


\d^ud‘^ 2(m-i) ||v|jC?°'M| 2(m-i) 

m— |q:| — 1 |q:| 


< ||n||m|VM 


^ m— |q:| — 1 


m—|ck| —1 


IH 


|c| 


\W\su\~-^\\Msu\\Z-\ 


< 


u 


3 

m' 


Using (10), we also estimate each term of the last snm of (12). For 0 < a < s, one has: 

(14) J d^u\v\s{d^ud‘^~°'u) < J \d^u\2\V{d°'ud^~°‘u)\2. 

We have |v(9"M(?^""n)|2 < \d°‘Vud^~°‘u \2 + |c?"mc?®""Vm| 2 - We estimate the hrst term exactly 
as previonsly. For the second term we obtain, again nsing Gagliardo-Nirenberg ineqnalities, 

< \d°‘u\ 2(m-i) |y~"Vn| 2(m-i) < ||n||m- 

|q:| —1 m— |q| 

We thns have obtained a Riccati-type differential ineqnality < G||m||^ which boils 

down to 

with a constant C independent of 5. This shows that the solntion can be extended to a time 
of existence T independent of 5 as well. Namely, we have the bonnd 


and so the critical time is T* 


\\u{t)\\m < 

(C'fnolU)-^. 


1 - Ct\\uo\\m' 


Step 4: Limit. For each 5 > 0, let be the solntion to (8) with the same initial data Uq. By 
the previons reasoning, us e C'([0,T); FT"*) nniformly in 5 for any hxed time T < T*. Let ns 
hx T <T*. Then, since is a Banach algebra, we estimate the right-hand side of (8) by: 

||ll(5| V|(5ri5 ~ IV |(5 ( 1 I 5 ) II m-1 ^ 111(51| m-1 ||ll(51| m || Us || m ~ II'^5 || m- 

This shows that us e C'^([0,T);iF™-'^) nniformly in S. 

We now tnrn to the convergence issne. Instead of relying on the Lions-Anbin compactness 
lemma, which only provides a limit for a snbseqnence, we show a more robnst convergence 
statement for the family us as 6 ^ 0. We claim that the family is a Canchy seqnence in 
G([0,T);L2). As a conseqnence of the interpolation ineqnality ||/||m' ^ ||/|lo~”*'^”^11/11 m 
for m' < m, it also means that the seqnence is a Canchy seqnence in any C{[0,T)] H^') for 
any m' < m. To prove onr claim, we need another estimate on the difference of operators 
|V| 5 . Let ns hx 5, £ > 0. The symbol of the difference is 

Sym(|v|,s - |V|.) = |e| /'"(e-^'^'^l - e-^^l«l)dA, 

DO 

which is bonnded nniformly in ^ by - s\. Therefore we obtain the following bonnd: 

(15) ||(|V|5-|V|,)n|U<C'|5-£|||n|U+2. 


Writing the eqnation for the difference, we obtain 

{Ue - Us)t = (w - Us)\y\eUe + M5(|V|e - |V|,5)ll£ + Usl^ls^Ue - Us) 

-iMe-Ms){ul)-Msiul-U^s)- 
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Testing with {u^ - us) we further obtain 

=y' {u6-U5y\y\eUe + J Us{Ue - Us){\V\e - \V\5)Ue + J Us{Ue - Us)\V\s{Ue - Us) 
- J (Ue-Us){\^\e-\V\s)iul) - J - M^)! V|5(m£ - M^), 

where in the last term we swapped onto (u^ -us). We see that the third term w^dvle - 
|V| 5 )n£(-Ue - Us) cancels with part of the last, and we have, using the same trick as above: 


J UeiUe - Us)MsiUe - Us) < -^ J U^\V\s{Ue - Us)'^ = -^ J (Ue - Us)‘^\V\siUe) < \u, - Us\ 


li, 


in view of the uniform bound on u^ in hf™. The rest of the terms can be estimated using 
(15). Note that m > (i/2 + 1 and being an integer, m > 2 for any dimension d > 1. So, 
Hm ^ Pg^ve 


A 

dt 


\Ue 


Us\l < C{\Ue - Us\l + - e\\Ue 




where C depends only on the initial conditions and other absolute dimensional quantities, 
but on e: and 5. Given that the solutions start with the same initial condition, the Gronwall 
inequality implies 

Wit) - usit)\2 < C\5 - e\ie^^ - 1) 

for all t <T. This proves the claim. 

So, the family us converges strongly to some u in all G([0,T); TT™'), m' < m. Moreover, 
dtUs converges distributionally to dtu, and in view of the uniform bound in it does so 

strongly in any This shows that the limit u solves (1) classically with initial condition 

uq. Passing also to a weak limit for a subsequence shows that u € G^([0, T);H"^) which is the 
space of weakly continuous H'^-valued functions. The argument to prove strong continuity 
in now follows line by line that of [20, Theorem 3.5] as we have all the same estimates. 
Finally, u e G^([0,T); follows as before for us, directly from the equation. 

Note that for the solution u that we constructed, the maximum principle proved earlier 
for Us still holds. The argument is the same, due to the positivity of the kernel. □ 


2.2. A Beale-Kato-Majda criterion. We now state the classical BKM criterion for our 
model. 

Theorem 2.2 (A Beale-Kato-Majda criterion). Suppose 

u e Ci[0,T)- H^iQ'^)) n C\[0,T)- H^-^Q'^)) 
is a positive solution to (1), where m > (i/2 + 1. Suppose also that 

(16) f \Vuit)\^dt < oo. 

Jo 

Then u can be extended beyond time T in the same regularity class. 

Remark 2.3. We will see that WluiWoodt < (X) is also a BKM criterion. 

The proof relies on a version of the log-Sobolev inequality a la Brezis [3] adapted to our 
setting. In order to state it, let us hrst recall some dehnitions. Let u ^ L"^. Then u admits a 
classical Littlewood-Paley decomposition u = E^-oo '^q- Let us denote the large-scale part by 
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u<o = 'ZqKo'^g and the small-scale part by m>o = Recall the classical (homogeneous) 

Besov norm 

|| m||rs =sup2®^|M„L 

-^r,oo ^ ' 


Lemma 2.4 (A log-Sobolev inequality). 

(17) |Vn|cx. + ||V|n|c« < \u\2 + I|w>o||b^^( 1 + log+ ||M||m) + 1- 

Proof. Bernstein’s inequalities imply |VMg|oo ~ ‘^^\uq\oo, uniformly for all q. In particular, we 
get 

oo 

|Vm|oo + ||V|m|oo ^ X! 2''|Mg|oo. 

q=—oo 

Then for g < 0 we use another Bernstein inequality: 2^\uq\^ < 2'^^^*^^\uq\2. Clearly, \uq \2 < \u \2 
and thus Eqi -00 2'?|Mg|oo ^ |n| 2 - For the small-scale component, we obtain, for every Q > 0, 


Q 


^2''KU = ^2^KU+ ^ 2'?KU<g||n>oiL 




g=0 


^Q\\u>oh 


q-Q+1 


+ X! 2'?F 

q^Q+l 


+^-m)29m 


Uqh 


Minimizing the above over Q (and recalling that m > d/2 + 1), the small-scale component is 
bounded by 


||n>o||si 1 + log 




Bi 


One observes next that -a: log a; < 1 on M+ and that logi/ < log+i/. Thus the small-scale 
component is overall bounded by 1 + ||m>o||bi (1 + log+ ||n||m), as claimed. □ 


We now turn to the proof of Theorem 2.2. 


Proof of Theorem 2.2. Let u be such that 

n +00 

(18) / ||m>o(OIIbi dt<+oo. 

According to the classical Bernstein inequalities, |VMg|tx, ~ 2'i\uq\oo, uniformly for all q, and 
by continuity of the Littlewood-Paley projections, |VMg|tx, ^ |Vn|tx). Hence (16) implies (18). 
Similarly, since there exist |v|-versions of Bernstein’s inequalities: ||v|nq|oo ~ 2i\uq\oo, the 
condition in Remark 2.3 also implies (18). 

Performing exactly the same estimates as on Step 3 of the proof of Theorem 2.1 but now 
with |v| instead of jv^l, we arrive at the following a priori bound 

M 

(19) dt\\u\l, ^ him E 

i=l 

where M depends on m and each of the /ij satishes 0 < /ij < 1. Indeed, in (12), the estimate 
of the symmetrized term of highest order gives /ij = 1. The estimate of the hrst sum gives 
Pi - l'^l/(^ “ 1) 0 < a < s and the specihc terms from the last sum are dealt with simply 

with /ij = 0. 

Combining (19) with the log-Sobolev inequality (17), we arrive at 

(20) dtWufm ^ I«|2||m|I™ + l|M>ollBi^_^(l + log+ ||m|U)||m||^ + ||m||^. 
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Applying Gronwall’s lemma twice, one gets a double-exponential estimate of the form: 

_log \\u{t)\\m _ 

log iwollm + t+ \u 0 \ 2 t + /q ||m>o||b^ 

Theorem 2.2 follows immediately. □ 




2.3. Prom local to global through regularity. We now study the question of the global 
existence and regularity of positive weak solutions starting from arbitrary L°°-data. 

Suppose that we are given a classical solution u e on the 

torus T'^ which is strictly positive m > 0. Let T* be its maximal time of existence. We will 
show that T* = 00 . Let us assume, on the contrary, that it is hnite. In what follows we apply 
the De Giorgi regularization result of [4] to our model. Since u is a classical solution, the 
formal passage from the u-equation (1) to the w = u^-equation (6) holds true. The active 
kernel k given by (7) is symmetric with respect to {x,y) and satishes 


( 21 ) 




< \k{t,x,y)\ < 


m 

\x - y^*^ 


for aW X i^y^ t> 0, and 

A(f) = Cdmax{\u{t)\^, |M‘Hf)|oo}- 

By the max/min principle, we see that A(t) can be replaced with A(0) = A, uniformly for 
all time, and thus depends only on |mo|oo- Equation (6) is exactly of the kind studied in [4]. 
It was natural in [4], in the context of an Euler-Lagrange problem, to assume the hniteness 
of the global energy, i.e. w € The main technical result of [4] however uses no such 

assumption and only requires w to have locally hnite energy (Gorollary 3.2 of [4] gives a global 
L°° bound in terms of the global norm, which in our case is not necessary as w remains 
bounded by the maximum principle). So, [4] applies verbatim to our periodic solutions, 
unfolded on the whole space M'^. Specihcally, the result states that there exists an a > 0 
which depends only on A and d such that for any 0 < to < E* we have w e x (tg, T*)) 

with the bound 

ll'?L'||c“f (T'^x(to,Tq) ^ ^*(^0, |wo|cx), A, d). 

Next, the Schauder estimates for integro-differential equations of type (6) recently obtained 
in [12] readily imply 


II^ II 1 4. ^ II^ II T* ' 1 U A? d) ^ C tg ; |^0 1 00 : A, d) . 

By the Beale-Kato-Majda criterion, tc, and hence u, can be extended smoothly beyond T*, 
resulting in a contradiction. 

A further application of the bootstrap argument of [12] readily implies higher regularity 
bounds (5) for w and hence for u. We thus have established the following result. 


Theorem 2.5 (Global well-posedness). Under the assumptions of Theorem 2.1, the solution 
exists globally in time. Furthermore, the solution is regularized instantly and satisfies the 
bounds (5). 
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2.4. Weak solutions from positive bounded data. The bounds (5) depend on the initial 
data only through the minimum and maximum value. This allows us to construct weak 
solutions from arbitrary positive data in L°° and for which similar regularization properties 
will hold. However it is not obvious how to restore the initial data and one should be cautious 
of the topology used for the limit t 0. 

Theorem 2.6 (Global weak solution). For any initial data uq e L°°(T^), uq > 0, there exists 
a global weak solution to (1) in the class 

X n n C'°(M+;L2). 

Its initial value uq is realized in the sense of the L°° weak* limit and in the strong sense. 
The energy is conserved, the momentum fjdu{x,t)dx is continuous on M+ and (4) is satisfied 
for any (t,t') € Furthermore, u satisfies the instant regularization estimates (5) and for 
all t > 0, the original eguation (1) is satisfied in the classical sense. 

Remark 2.7. If uniqueness was to fail in Theorem 2.6, it could only do so at t = 0. However, 
the continuity of the momentum at t = 0 prevents any concentration of the H^l‘^ norm in our 
weak solution. 

Proof of Theorem 2.6. Let Uq e L°° be such a positive initial condition. We start the con¬ 
struction by taking standard mollihcations (no)e of '^o- In view of Theorem 2.5, there exists 
a global classical solution from each of those mollihed initial conditions. Let us £x T > 0. 
In accordance with (5), the family {ufi^ satishes uniform regularity starting from any F > 0. 
Invoquing the Arzela-Ascoli compactness theorem together with Cantor’s diagonal argument 
to successively reduce the value of to, we can pass to the limit and hnd a classical solution 
on the interval (0,T]. In addition, as we already pointed out in the introduction, the en¬ 
ergy equality on the whole time interval [0,T] combines with the momentum law to ensure 
that «£ belongs to L°°([0,T) x n L2{\D,T)]H^I‘^) uniformly. As L°° n is an algebra, 
= u‘1 enjoys the same uniform property. Passing to the weak limit in the Sobolev space, 
we conclude that 

(22) u,w ^ L°°{[0,T) n L\[0,Ty,H^I^). 

For all t > 0, since we also have a limit in the classical sense, w{t) = u{t)'^. The only 
remaining problem is to restore the initial data and a weak formulation of the equation near 
t = 0. Indeed, once a solution is constructed on [0,T], one can invoque Theorem 2.5 one 
last time, starting with the smooth initial data u{T) and therefore claim the existence of a 
global solution on M+. 

Let us restore the initial data for w hrst. Let us write (6) in the weak form, for the 
smoothed solutions: for any (j) e C'°°([0,T) x T'^), 

(23) / wJx,t)(f){x,t)dx- / wJx,t))4){x,t))dx- / / We(x,t)dt(j)(x, s)dxds 

JTd Jjd Jo JTd 

= ^ f f Ks^^^y){we{y,s)-We{x,s)){(t){x,s)-(t){y,s))dxdyds. 

2 Jo Jt'^ 

Let us write (23) as A,, - with the respective designation of each term. As 

we pass to the limit £ 0, we trivially have -> Aq. The convergence Bq holds 

because (uo)e ->■ Uq strongly in and hence Ws{t = 0) = (uq)^ ->■ Wq in and thus weakly. 
The convergence Co results from the uniform L°° bound on We (in an arbitrary small 
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neighborhood of t = 0) joined with the uniform convergence on the rest of the time interval. 
Finally, one splits the last term as: 

(24) r+ 

Jo JO 

We have D'^^hy classical convergence and, by Cauchy-Schwarz and (22), 

(25) \^e,5\ - II'*^£|Il 2([0,(S);R1/2) ||0||£2([o,5);iyi/2) ^ C\f5 

uniformly in £. This shows that Dg, Dq. We have thus proved that (23) is satished in the 
limit as e 0, i.e. for w itself. Taking then 0 independent of t and passing to f 0 shows 
that w{t) Wq weakly* in L°°. 


At this point let us make a cautionary remark because u{ty does not imply u(t) uq 
in general. A simple example is provided by the sequence «„ = 1 + where rn{x) = 
sign(sin2”7rx) are Rademacher functions. Then > 0, 1 and yet -^ | 12 _ A 

progressively mollihed sequence (Mn)i/n" would provide a similar counter-example in the 
class C°°, as in our case. However, one can observe on this example that the function y^S/d, 
whose square is the limit of squares, dominates the limit of itself, which is 1. This is true 
in general. 


Lemma 2.8. Suppose that a sequence of functions {«„} c bounded away from zero, 
enjoys both limits Un u' and ^ in the weak* topology. Then uq>u'. 


Proof. The proof of this lemma is simple. Let 0 > 0 be a test function. Then trivially 

J {Un-UoYcfdx > 0, 

for all n. Let us expand, 

J {un-Uof(j)dx = J (ul-2unUo + ul)(j)^ J {ul - 2u' Uq + ul)(j) = 2 J Uo{uo-u')(j)>0. 
Since Mq > 0 and the above holds for an arbitrary 0 > 0 in L^, the lemma follows. 


□ 


Let us go back to restoring the intial condition for u. Reverting to the original equation (3), 
its weak formulation for the smooth sequence reads, with test functions independent of t: 


(26) / Ue{x,t)4){x)dx - / Us{x,0)(p{x)dx = 

JTd- 


(j)(x){u,;(y,s)-Ue{x,s)yKper{y-x)dydxds 

2 Jo Jt^ 


1 rt 
+ — 

2 


f [ Ue{x,s){(t){x)-(t){y)){ue{y,s)-Ue{x,s))Kp^,{y-x)dydxds. 

Jo Jjd- 


Passing to the limit on the left-hand side and in the last integral on the right presents no 
difficulty as one can use estimates similar to (24)-(25). 

However there are no a-priori bounds that guarantee the smallness near the time t = 0 of 
the hrst integral on the right-hand side. Specihcally, a possible concentration of the 
norm near t = 0 (or equivalently by (4), an initial discontinuity in the hrst momentum) could 
prevent (26) from withstanding the limit. 
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Using only the positivity of the first integrand of the right-hand side, we obtain instead, 
for all 0 > 0: 


(27) / u{x,t)4){x)dx - / Uo{x)(j){x)dx > 

JT'i JT'i 

u{x,s){(j){x) - (j){y)){u{y,s) 


1 r. 

2 Jo 


u{x, s))K-pej.{y - x)dydxds. 


In the limit t ->■ 0, the right-hand side of the previous inequality vanishes by (25), and hence, 
any weak*" limit of a subsequence of (M(t))t>o would converge to a function u' satisfying 
u' > Mo a.e. On the other hand, by the lemma above, uq > u'. This proves that the weak* 
limit u{t) Mo holds as t -> 0 for any subsequence and therefore also as a function of the 
continuous time parameter. In particular, testing this weak* limit with 0 = 1 ensures that 
the momentum Jjau{x,t)dx is continuous at t = 0. 


Let us now recall that the norm is also continuous at time zero and therefore overall 
preserved. Indeed, one can for example use 0 = 1 as a test function in the weak* limit 
w(t) Wo or directly notice that \ue(t)\2 = |(mo)£|2 for e > 0 and that this identity passes 
to the limit c ^ 0, on the left-hand side because of the uniform convergence at time t > 0 
and on the right-hand side because it is a standard property of mollihed sequences. Next, 
as Mo e c is an admissible test function for the weak* convergence u(t) Mq, one gets 

|M(t) - m( 0)|2 = 2|moP - 2 f u{x,t)uo{x)dx ^ 0. 

Therefore m e C'(M+;L^). 


Finally, to restore the weak formulation of the original equation, we hrst write it on a time 
interval [toT] with to > 0, where it is also satished in the classical sense: 


(28) / u{x,t)(j)(x,t)dx - / u{x,to)4>{x,to)dx - / / u{x,s)dt(j){x,s)dxds = 

Jt^ Jto 

(j){x, s)(u{y, s) - u(x, s))^Kper(y - x)dydxds 

f f u{x,s)((j){x,s)-(j){y,s)){u{y,s)-u{x,s))Kper{.y-x)dydxds. 
Jto MT'* 


2 Jto MT'i 
1 

+ — 

2 Jto 


Passing to the limit as to 0 is now within reach. On the left-hand side, the continuity 
ensures the convergence of the middle term while the L°° bound tackles the third term. On 
the right-hand side, the last integral is controled in the same fashion as (25) because of (22). 
The troublesome hrst term on the right-hand side is bounded effortlessly; for any 0 < U < fo^ 


(29) 


f f (t){x,s){u{y,s)-u{x,s)yKper{y-x)dydxds <C f 

Jti Jt^ Jti 


Ms)f^y^ds. 


Applying (4) on [fi,fo], where the equation is already satished in a classical sense, the right- 
hand side equals C fjdU{x,to) - u{x,ti)dx. As the momentum is continuous, the right-hand 
side is arbitrarily small as to ^ 0, uniformly for ti e [0,to]- One can thus hrst let ti 0 in 

(29) and then pass to the limit to 0 in (28). Finally, taking 0 = 1 in the weak formulation 
of the equation shows that (4) still holds at t = 0. This hnishes the complete construction of 
a weak solution on [0, oo) and the proof of Theorem 2.6. □ 
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2.5. Blowup in finite time for some smooth negative data. In view of the time 
reversibily property (TR) mentioned page 3, if m is a positive solution to (1), then -u{t* -t) 
is a negative solution. Thus starting with positive data Uq e L°°(T'^)\C(T'^) we obtain a 
solution u which becomes C°° instantaneously. Then -u{t*) serves as negative initial data 
that develops a singularity at time t = t*. 

Corollary 2.9 (Finite time blow-up). For any t* > 0, there exists a negative initial con¬ 
dition uq e C'°°(T'^), Mo < 0 and there exists a classical solution to (1) that develops into a 
discontinuous solution at time t* i.e. u{t*) € L°°(T'^)\C(T'^). 


3. Long-time behavior 


In this section we will show that the long-time dynamics of the model converges to a 
constant state consistent with the conservation of energy, namely, 


u(t, x) 


t^ + oo 


II'»oIIl2 

vTS^' 


As the solution is squeezed by the maximum and minimum principles, it is expected that 
such behavior would be a consequence of a steady decay of the amplitude. We will show 
indeed in Lemma 3.1 that the amplitude tends to zero exponentially fast. However, this 
first property alone still leaves room for non-trivial oscillations like a non-vanishing |Vm|oo- 
We will exclude such residual oscillations using the method that Constantin and Vicol [7] 
developed recently to prove, among other things, the global regularity for the critical SQG 
equation. In short, the long-time dynamics of our model is a convergence to a constant state, 
both at large and small scales. 

Subsequently, we use the following notations: 


m(t) = minM(a:, t), M(t) = maxu{x,t), A{t) = M{t) - m{t) 
and mo = m(0), Mq = M(0), Aq = A(0). 

Lemma 3.1 (Exponential decay of space oscillations). Let u € L°°(T‘^ x [0, cx))) be a weak 
solution to (2) with a positive initial condition as stipulated in Theorem 2.6. Then A(t) < 
holds for all t > 0 with some c> 0 independent of u. 


Proof. This proof relies on an idea from [11]. Let x{t) denote a point where M(t) = u{x{t),t) 
and let x(t) denote a point where m(t) = u{x{t),t)- Let us unfold the 27r-periodic solution 
on ; without loss of generality, we can assume that x,x^ [-7r,7r]'^. The gradient \/u{-,t) 
vanishes both at x{t) and x{t). Then we have, for instance in the viscosity sense of Crandall- 
Lions, 


M'{t)< J u{y){u{y)-u{x))K{y-x)dy 

<mo / {u{y)-u{x))K{y-x)dy 

J\y-x\>l,\y-x\>l 

< mo / {u{y) - u{x))m\n{K{y - x),K{y - x)) dy. 


We dropped the reference to time in the right-hand sides for readibility. Similarly, we have 
m'(t) > mo / {u(y) - u{x)) min (K(y - x), K(y - x)) dy. 

J\y-x\>l,\y-x\>l 
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This implies that 

A'it) <-moA{t) I min {K{y - x),K{y - x))dy <-cmoA{t) 

where c = yiy|>i+^^ follows by Gronwall’s lemma. □ 

Theorem 3.2 (Asymptotic behavior at all scales). Let u € x [0, oo)) he a weak 

solution to (1) with a positive initial condition as stipulated in Theorem 2.6. Then m{t) is 
increasing, M{t) is decreasing and both A(t) and |VM(t)|oo decay to zero exponentially fast. 
Conseguently, the solution converges to a constant at an exponential rate. 

Remark 3.3. Combining this asymptotic with (4) allows us to compute the total increase of 
momentum as the defect in the Cauchy-Schwartz inequality for uq : 

- J^^uo(x)dx. 

Lemma 3.1 indicates that the stabilization of the largest scales starts right away. The 
following proof will establish that the time-scale of the transitory regime that precedes the 
stabilization of the lowest scales does not exceed T* 2 ; ^ log+(Ao/Cmo). 

Proof. Let us take the gradient of the integral form (2) and multiply by Vm. After elementary 
manipulations we obtain pointwise (the integrals being understood as principal values): 

^dt\\/u{x,t)\‘^ + \/u{x) ■ J{yu{x)-Vu{x + z))u{x + z)K{z)dz 
= Vu{x) ■ J{u{x + z) - u{x))Vu{x + z)K(z)dz. 

Inside the hrst integral let us use the identity 

Vm(x) • (Vm(x) - \/u{x + z))u{x + z) = -[|VM(a:)P - |VM(a: + z)\f^u{x + z) 

(30) I 

+ -|Vm(x) - Vm(x + z)f‘u{x + z). 

Assuming that x is the point of the maximum of IVmI, the first term on the right-hand side 
of this last identity is non-negative, while in the second term we can simply use the bound 
from below u{x + z) >m. We thus obtain the estimate 

dt\Vu{x,t)\^ + mo y" \Vu{x + z) - Vu{x)\^K{z)dz 

<2 / SJu{x) ■S/u{x +z){u{x +z)-u{x))K{z)dz. 


Let us denote the integrals respectively by / and J: 

(31) dt\Vu{x,t)\^ + < J. 
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We split the integral J depending on whether \z\ < p or | 2 :| > p and write J = J<p + J>p. We 
estimate J>p after rewriting it in the following form 


J>p ~ 


2\Ju{x)- / Vz{u{x + z) - u{x)){u{x + z) - u{x))K{z)dz 

J\z\>p 

\/u{x) ■ / 'Vz{u{x + z) - u(x)yK{z)dz 

J\z\>p 


{d+l)Vu{x)- f z{u{x + z) - u{x)Y + Vu{x) ■ f Vz{u{x + z) - u{x)) 

J\z\>p \z\^'^'^ J\z\=p 


,Cdda{z) 


P' 


d+l 


Thus, 


p2 ■ 


\J>p\ < C2\Vu[ 

As to J<p, we use the first order Taylor formula for the increment of u\ 

J<p = 2 f p.v. ( r S/uix) ■ Vu(x + z)S/u(x + Xz) ■ zK(z)dz\d\ 

Jo \J\z\<p } 

= 2 f f SJu{x) ■S/u{x + z){SJu{x + Xz) - S/u{x)) ■ zK{z)dzdX 
Jo J\z\<p 

+ 2 f f S/u{x) ■ {SJu{x + z)-SJu{x))SJu{x) ■ zK{z)dzdX 
Jo J\z\<p 


= + j2 

^ <p ^ ^ <p- 


We get an upper bound of J^p by Cauchy-Schwarz: 


\J\p\<2\Vu\l J\ 


|Vm(x + z) - Vm(x)| 


i\z\<\p 


Cddzj dX 

\/Cd 


2|Vm|L [ |VM(a; + z) - Vu{x)\\/K{z) - -^dz 

J\A<p \z\ 2 


< 2\Vu\'i,y/cdUJdIp < p. 

4 rriQ 

The estimate for is completely analogous (without the A-integral). Thus, 


T ^ "^0 r ^ 
J<p < —^I + C 4 


Wu\t 

mo 


-P- 


Incorporating the estimates already obtained back into (31), we arrive at 

/ M9 LUn ^ IVmII 

dt\Wu{x,t)\^ + —I < C2|Vm|oo^ + C 4 — —p. 




mo 


The choice p= {moA'^y/^l\S/u\oo optimizes the right-hand side: 


( 32 ) 


dt\X/u{x,t)\^ + ^I<C5l — ] |vm|: 

2 Vmo/ 
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Next, we estimate I from below using an argument analogous to that of P. Constantin 
and V. Vicol [7] but in which we incorporate the amplitude. For an arbitrary r > 0, we write 


1 = 


[ \\/u{x) - S/u(x + z)\‘^K(z)dz > [ \S/u{x) -\/u{x + z)\‘^K{z)dz 

J J\z\>r 

I \Vu{x)\‘^K{z)dz - 2Vu{x) ■ f Vu{x + z)K(z)dz 

J\z\>r J\z\>r 


therefore 


I> 

I >c& 
= C6 


P|: 

|Vm(x)P 

r 

|Vm(x)P 


^ / X r Vziu(x + z) - u(x)) 

- 2Wu{x) ■ J ^ ^ ^U ,dz 


'\z\>r 


-2{d+l)\/u{x) ■ f z(u{x + z) - u(x)) 

J\z\>r \z\ 

2\Juix)- f i'Ju(x + z) ■ 

J\z\=r 


It follows 


|VM(a;)p . .A 

I > Cq- - - - Cy\\/u{x)\ — 


and choosing r = provides 


(33) 


/ > 2ci 


|Vn(x,f)|3 
‘ A 


dtWu\l + Cs'^\Vu\l<cJ — Y |V 
A \mo/ 


As a: is a point of maximum of \/u{t), combining (32) with (33) yields 
(34) 

In view of Lemma 3.1, there exists a time T* ^ ^ log^(Ao/C'mo) such that for t>T*, 

(Ay mo 

C5 ^^8 —■ 

Vmo/ 2A 


This implies that for t>T* 
or equivalently. 


5i|Vn|L + C8^|Vn|L <0 


5t|Vn|oc + C8^|Vn|L <0. 


Using the precise estimate from Lemma 3.1, we further get 

dt\Vu\oo + < 0 . 

4/lo 


This hnally implies for t>T* 
|Vn(f)|oo < 


|Vn(T*)|c 


|Vn(T")Ue-^™o* 


C9|V4t(T*)|oo /^cmnt _ ^cmoT*\ _|_ C9|Vn(T*)|oo /_ ^-cmo{t-T*)\ 

Aq t ' Aq ' ' 

The proof of Theorem 3.2 is now complete. 


□ 
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4. Local well-posedness in analytic classes 

In the last section we plan to explore numerically what happens, in terms of existence or 
blow-up, to solutions when the initial condition is unsigned or even negative. Since for such 
data we do not even have a local existence result in Sobolev spaces, we hnd it necessary to 
prove a generic local existence result, at least in analytic classes. Since our numerical data 
obviously has a finite Fourier spectrum, it will ensure the minimal solid ground required to 
run the numerical scheme. 

We rely on a not so well-known hxed point theorem by Nishida and Nirenberg, [23]. Let 
us recall it. Let us define the following spaces, for p > 0, 

(35) Xp = {u-.\ul= ^ e\^\P\u{k)\<oo), 


(36) 




Notice that Xp is a Banach algebra. Also notice that for any p' > p one has 


(37) 

for some absolute (7 > 0. 


C 


|Vu|p = ||V|m|p < -|m|p'. 


p - p 


Theorem 4.1 (Nishida and Nirenberg, [23]). Suppose a functional F = F{u) satisfies the 
following condition: F -.Ypi Yp for all p' > p, and for any R> 0 there is a Cr such that for 
all \\u,Ui,U 2 \\p' < R and for all p < p' the following holds: 

(i) IIUUls 

(ii) ||F(!1 i) -f(!l2)lli ^ll“l 
Then the Cauchy problem 


lu(0) = 


is locally well-posed in the following sense: for any Uo e Yp^ and any p < p', there exists T 
and a unigue solution u e C^{{-T,T)]Yp) to (38). 

In our case, F{u) = u\y\u - |v|(m^). Let us fix an i? > 0 and assume that || m || p ' < R and 
p < p'. Using the algebra property of Yp and (37) we estimate (absolute constants may 
change from line to line) 


and 


\ F ( u)\p < |m|p||v|m|p + 


C 

P'-P 



CR , , CR 

- M o' < - 

P'-p P'-p 


U 


p'5 


\YF{u)\p < |v(m|V|m)|p + 2||v|(mVu)|p 

C C 

-(|m|p'||V|m|p' + |m|p'|Vu|p') < 


p'-p 


p'-p 


Wfp < 


CR 

p'-p 


M\p'- 


Thus, ||F(m)||p < ^^||m||p'. The Lischitzness condition is verified analogously. 
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5. Numerical simulations and perspectives 


In this last section, let us illustrate briefly with some numerical simulations the main 
analytical results of our paper, in dimension d = 1. 


5.1. Smooth positive data. Positive smooth data satisfy the assumptions of our main 
Theorem. We computed the solution using a direct integration of the singular-integral form 
of (2). Given a time step r > 0 and a spatial discretization sequence {xi)i<i<N of 2 ; [-vr, vr] 
with a uniform mesh size 6 = Xj+i-Xj, the value Ui^k representing u{kT,Xi) is computed with 
a forward numerical scheme: 


(39) 


Uk+l,i — Uk^i + t5 

j*i 


{UkJ Uk,i)Uk,j 
4sin2(^) 


The code runs well as long as one makes sure to respect a CFL (r is small enough) and that 
each term stays well within the limit of exact computer arithmetic (to compute accurately 
the compensations i.e. the principal value). 

Figure 1 shows the evolution from uo{x) = 2 + sin(x) + ^cos(5a;) and the profile of the 
solution for various time-slices. As expected in view of Theorems 2.5 and 3.2, the solution 
converges rapidly towards the theoretical limit 


WuqWl^ 



2.1319. 


The maximum and minimum principle are clearly visible. 




Figure 1. Evolution from various data. Left: mo( 3^) = 2+sin(a;) + ^ cos(5a;). 
Right: uo{x) = 1/2 if a; e [-vr, 0) and 1 - ^ if x e (0, vr]. 


5.2. Evolution from discontinuous data. Figure 1 shows also the evolution from a pos¬ 
itive initial condition with a single discontinuity: 


uo{x) 


[ 1/2 if X e [-vr, 0) 
jl-l^ ifxe(0,7r]. 


Even though the numerical scheme (39) is not particularly refined, we did not observe any 
numerical instability. This observation is consistent with the instant smoothing effect of our 
Theorem 2.6 and is a timid but positive point to support the conjecture that only one weak 
solution arises from positive bounded data. We pushed the test further and Figure 2 shows 
the evolution from non-smooth numerical data cooked up to simultaneously be non-derivable 
at ±71 and -1, have a jump at 0 followed by a chirp-like oscillation ^/xsin(l/x). The energy 
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spectra (modulus of Fourier modes, in log-log scale) of various time-slices attest to the reality 
of the smoothing effect over 2 decades and beyond. 



Figure 2. Evolution from non-smooth data (non-derivable at ±7r and -1, a 
jump at 0 followed by a chirp Aya:sin(l/x)) and energy spectrum. 

5.3. Evolution from negative data. Figure 3 shows the evolution starting from 

3 

uq{x) = -2 - sin(x) - — cos(5x), 

i.e. the opposite of the positive data of paragraph 5.1. The first time-slices remain smooth 
as one could expect from the local well-posedness in analytic classes of Section 4. The 
(AMP) principle is clearly visible as well. The last time-slice shows the emergence of the 
first singularity: a highly oscillatory singularity forms around the minimum value of Uq, 
while the rest of the solution remains relatively stable. A secondary instability is also barely 
discernible near the second lowest value of the data. 























































































































































22 


CYRIL IMBERT, ROMAN SHVYDKOY, AND FRANCOIS VIGNERON 




Figure 3. Evolution from negative data and evolution of the Fourier Spec¬ 
trum showing the emergence of a singularity. 
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5.4. Evolution from unsigned data with a positive average. Let us write the Non¬ 
local Burgers equation (1) in a form that separates the average value of the solution from 
its fluctuations around it. To this end, we introduce the zero-average function 


v(t, x) 



The average value can be recovered from v by the momentum law (ML) 

f u{t,y)dy= f uo{x)dx+ f ||n(r)f^i/2dr. 

JTd JTd Jo ^ 

Note that the average of u is globally bounded from above and below by: 

7—77 f uo(x)dx<—^ f u{t,y)dy f uoixYdx 

\Td\JTd > IT'^I V ^ ’ 

The zero-average fluctuations v evolve according to 



(40) 



Therefore, if the initial condition has a positive average, equation (40) contains a uniform 
smoothing term on the left-hand side, while the right-hand side is a non-local, quadratic but 
energy preserving non-linearity of order 1 i.e., formally: 

Ik(0lli2 + 2 ^ J^^u{T,y)dy'^ \\v{r)\\fjy2dT = \\vo\\l2. 


Figures 4, 5 and 6 show numerical simulations starting from different unsigned data but 
with a positive average. For the first one, the slight incursion in negative values does not 
seem to affect the solution. The second one shows a localized transient singularity (which 
might have been smoothed out by the finite size of the mesh). The Fourier spectra at various 
time slices show the emergence and the subsequent vanishing of the singularity. 

Figure 6 displays an interesting feature of those transient instabilities: they appear to be 
localized along the negative local minima of the data. Their sequence of development seems 
also to start from the one with the highest absolute value and then proceed in order towards 
the ones of lower magnitude. 


Remark. In its form (40), the Non-local Burgers equation reveals some similarity with the 
incompressible Navier-Stokes equation that can be written as 

dtu - uAu = -u ■ Vn - Vp 

with divM = 0. We already pointed out the similarity of our model with the Euler part 
of NSE, especially if one computes the pressure in term of u. For the fluctuation v we 
also consider that the active “viscosity” term on the left-hand side of (40) is similar to the 
kinematic viscosity in NSE. Despite this similarity with the notorious counterpart, we were 
able to prove global existence for our model. Of course, in our case this boils down to knowing 
extra structure, e.g. the maximum principle. This once again highlights the importance of 
finding additional structural properties of the NSE in any attempt to approach its celebrated 
global regularity problem. 
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Figure 4. Evolution from mq = 1 + sin(a;) + cos(20a;). 
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Figure 6. Evolution from Uo = ^ + \ sin(19x) + | cos(20a:). 

5.5. Global simulation with a finite element method. To minimize the effects of time- 
discretization, we did numerical simulations of (2) using finite elements, based on piecewise 
linear functions. The resulting problem is an ODE of the form AU' = J{U,U) where U{t) e 
M"’ represents the value of {u{t,Xi))i<i<n, ^ is a Toepliz matrix and J is an n-dimensional 
vector of quadratic forms. As energy is also preserved for the ODE {^UAU = Cte), it is 
globally well posed for all times f e M (forwards and backwards). 

As a happy coincidence due to the special structure of the non-linearity, diagonalizing the 
Toepliz matrix A drastically simplihes each quadratic form in J, which then allows for a 
direct computation of the solution. 

Figure 7 shows two solutions obtained in that manner with 20 grid-points. The interesting 
point of that simulation is that the flow seems to connect a negative constant ground state 
to a positive one and that all the oscillations occur in the non sign-dehnite region. The 
level-sets of that region are shown above the 3D perspective. 

5.6. Perspective: the Frozen Non-local Burgers model. As the previous numeri¬ 
cal simulations have shown, the most interesting and mathematically challenging region 
is when u is non sign-dehnite. We feel that it is worth mentioning that the Non-local Burg¬ 
ers equation (1) comes with a natural twin in which the average value is frozen instead of 
memorizing the past norm. The Frozen Non-local Burgers equations reads: 

(41) = 

This equation should be compared with the previous equation (40) of the huctuations v of 
the solution of (2). Formally, (41) conserves momentum fjav(t,x)dx = Cte. The previous 
intuition suggests that it would be most interesting to study it for vanishing average data, 
in which case it formally conserves the L^-norm too. 
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Figure 7. Global FEM simulations. 
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